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• In this paper we consider the Extra-solar Planetary Systems recently discovered in our 

^1^' Galaxy as potential sources of gravitational waves. We estimate the frequency and 

characteristic amplitude of the radiation they emit due to the orbital motion, using 
the quadrupole formalism. In addition, we check whether the conditions needed for the 
resonant excitation of the /- and g-modes of the central star can be fulfilled. By a Roche- 
ts ■ lobe analysis, we show that there could exist systems in which the low-order g-modes 

could be excited, although this docs not happen in the systems discovered up to now. 



o 



■ 1. Introduction 

' Since 1992, when the first Extra-solar Planetary System (EPS) was discoveredB, 

a number of such systems, composed of one oil more planets orbiting around a 
main sequence solar-type star has been observedB, and it is conceivable that many 
others will be found in the near future. Some of the discovered EPS's are well 
characterized, since the mass and radius of the main star, the mass of the observed 
planets and their orbital parameters can be deduced from observations. With this 
information it is possible to attempt a first estimate of the characteristics of the 
gravitational radiation emitted by these systems, that, being at a distance of a few 
tens of parsecs, are very close to us. We shall select a set of EPS's for which all 
needed information is available, and compute the gravitational signal emitted by 
each couple star-planet due to the orbital motion, using the quadrupole formalism. 
The estimated frequency and amplitude of the waves impinging on Earth will be 
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compared with the emission of the binary pulsar PSR 1913+16Q'Q'cl. 

We shaU then consider a further mechanism of gravitational emission which may 
play an interesting role in these systems, i.e. the resonant excitation of the modes 
of oscillation of the star induced by tidal interaction with the orbiting planet. This 
problem has been studied in the literature for solar type stars with a planet, and 
for neutron star-neutron star close binary systems by using different approaches, 
based essentially on the deformation induced by the dynamical tides raised by the 
companionEl'Q^Q. The emission of gravitational waves associated with the resonant 
excitation of the modes of a star has been studied by Kojimafl. He considered a 
small mass in circular orbit around a compact star, and integrated the equations de- 
scribing the perturbation of the star excited by the small mass, in general relativity. 
He showed that a sharp resonance occurs if the frequency of the fundamental mode 
of the star is twice the orbital frequency, and that the characteristic wave amplitude 
emitted at that frequency can be up to 100 times larger than that evaluated by the 
quadrupole formula. In this view, it is interesting to check whether the conditions 
of resonant excitation can be fulfilled in planetary systems for some modes of the 
central star, and in particular if this is the case for any of the EPS's listed in table 

1. This will be done in section 4. In addition, we shall discuss the possibility for 
a planet to get sufficiently close to a star to excite the lowest order g-modes, and 
possibly the fundamental one, without being disrupted by tidal forces. We will show 
that for solar type stars the resonant excitation of the low-order i^-modes may, in 
principle, be possible. 

2. Main characteristics of the extrasolar planetary systems 

Among the EPS's discovered up to now, we have selected those for which the 
parameters of the central star and of the planets, which we need to estimate the 
gravitational emission, have been determined with sufficient accuracy. These pa- 
rameters are tabulated in table 1 and 2. Here and in the following, data will be 
given with the corresponding errors, when available in the literature. In the first 
column of table 1 we list the selected EPS's with the corresponding bibliography, 
and in column 2 the spectral class of the central stars. Most of them belong to a 
class similar to that of the Sun, which is a G2 V star. 

In column 3, 4 and 5 we tabulate, respectively, the mass of the centraLstar, A'/*, 
its distance from Earth, D, - taken from the Extrasolar Planets Cataloga- and the 
quantity which will be used in section 4. It should be mentioned that the 

presence of planets has been discovered mainly by using accelerometric techniques, 
that is to say, by measuring the variations of the radial velocity of the star caused 
by its gravitational interaction with the planet. These techniques do not allow to 
determine the mass of the planet, Mp, but only the product AfpSini, where i is 
the angle between the line of sight and the normal to the orbital plane. In column 6 
we list the values Mp sin i for each planet, normalized to the mass of the central star. 
From the data of table 1 we see that the closest system is at 4.70 pc, the farthest 
at 45 pc, whereas the mass of the central star ranges within [0.3, 1.37] Mq. 
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Table 1 

Earth, D, the factor 



The spectral class of the central star, its mass, M*, its distance from 



(rad-s 



and the ratio between the lower limit of the 



estimated mass of the planet, MpSini, and M^,, are tabulated for a selected set of 
EPS's. Strictly speaking, the bodies orbiting the stars listed below HD 114762 are 
classified as brown dwarfs rather than as planets. 



Star 




Spectral CI. M^{Mq) D (pc) 



IQ^ MpSini/M^ 



HD 75 
51 PeptH 
V Andl"! 



55 CnpM 
p CtBB 
HD 2102'ny 
16 CytBEl 
Gl 8761 
47 U: 

14 Hi 

Gl 

r Bool 
HD lfiM4 
70 ViiO 
HD 114762 
HD 110833!^ 
HD 11275«l" 
HD 2958'^^ 
HD 28375, 
HD 89707t4 
HD 21758011 




GO V 


1 


05 


28.94 


5.4 


3 


8 


10 


-4 


G5 V 


1 


05 ±0.09 


15.36 


5.1 ±0.7 


4 


2 


10 


-4 


F8 V 


1 


34 ±0.12 


13.47 


3.7 ±0.5 


5 


3 


10 


~4 












1 


5 


10 


~3 












3 


3 


10 


-3 


G8 V 





95 ±0.10 


12.53 


6.6 ± 1.0 


8 


8 


10 


-4 


G2 V 





89 ±0.05 


17.43 


3.8 ±0.5 


1 


1 


10 


-3 


GO V 





92 


21.29 


5.0 


1 


3 


10 


-3 


G5 V 





96 ±0.05 


21.62 


5.0 ±0.6 


2 





10 


~3 


M4 V 





3 


4.70 


13 


6 


7 


10 


~3 


GO V 


1 


01 ±0.05 


14.08 


4.5 ±0.5 


2 


2 


10 


-3 


KO V 





8 


18.15 


8.0 


3 


9 


10 


-3 


Kl V 





79 


10.91 


8.0 


4 


4 


10 


-3 


F7 V 


1 


37 ±0.09 


15.60 


4.7 ±0.6 


3 


3 


10 


-3 


G5 V 





94 


37.88 


7.2 


5 


1 


10 


-3 


G5 V 


1 


01 ±0.05 


18.11 


2.5 ±0.2 


6 


9 


10 


~3 


F9 V 





75 ±0.15 


40.57 


4.0 ± 1.0 


1 


3 


10 


-2 


K3 V 





75 


17 


8.2 


2 


2 


10 


-2 


KO V 





8 


16.5 


8.0 


4 


2 


10 


-2 


G2 V 


1 





45 


6.1 


3 


8 


10 


-2 


K2 V 





75 


16.5 


8.4 


6 


4 


10 


-2 


Gl V 


1 


2 


25 


6.1 


5 





10 


-2 


K4 V 





7 


18 


9.1 


8 


2 


10 


-2 
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Table 2: Orbital parameters of the planets orbiting around the stars listed in table 
1 (1 AU = 1.495978706-10^3 cm). 



Star 


P 


a (AU) 


e 


HD 75289 


3.5103 ± 0.023 


0.046 


0.054 ±0.026 


51 Peg 


4.2293 ±0.0011 d 


0.0520 ±0.0015 


0.012 ±0.01 


V And 


4.6170 ± 0.0003 d 


0.598 ±0.018 


0.109 ±0.04 




241.2 ± 1.1 d 


0.83 


0.18 ±0.11 




1266.6 ± 30 d 


2.50 


0.41 ±0.11 


55 Cnc 


14.648 ± 0.0009 d 


0.1152 ±0.0038 


0.051 ±0.13 


p CrB 


39.645 ± 0.088 d 


0.2192 ±0.0042 


0.15 ±0.03 


HD 210277 


437 ±25d 


1.097 


0.45 ±0.08 


16 Cyg B 


804 ± 11.7 d 


1.66 ±0.05 


0.634 ±0.082 


Gl 876 


60.85 ±0.15 d 


0.21 ±0.01 


0.27 ±0.03 


47 Uma 


3.0 2/ 


2.08 ± 0.06 


0.03 ±0.06 


14 Her 


1619 ± 70 d 


2.5 


0.3537 ±0.088 


Gl 86 


15.83 d 


0.11 


0.05 


T Boo 


3.3128 ± 0.0002 d 


0.0413 ±0.001 


0.018 ±0.016 


HD 168443 


57.9 d 


0.277 


0.54 


70 Vir 


116.6 ±0.01 d 


0.47 ±0.01 


0.4 ±0.01 


HD 114762 


85.03 ± 0.08; d 


0.34 ± 0.04 


0.35 ±0.05 


HD 110833 


270.04 d 


0.8 


0.69 


HD 112758 


103.22 d 


0.35 


0.16 


HD 29587 


3.17 y 


2.5 





HD 283750 


1.79 d 


0.04 


0.02 


HD 89707 


298.48 d 


0.91 


0.93 


HD 217580 


456.44 d 


1 


0.52 



In table 2 we tabulate the orbital parameters of the planets orbiting around 
the stars listed in table 1, i.e. the orbital period P, the semimajor axis a, and the 
eccentricity e. From the data in column 2 we see that in some cases the period 
is very short, indicating that the planet gets very close to the central star. For 
instance, it can be as short as 1.79 days for the planet orbiting around HD 283750. 

3. Gravitational wave emission due to the orbital motion 

We shall first evaluate the amount of gravitational radiation that each couple 
star-planet emits because of the orbital motion. Since we are not interested in 
the detailed form of the signal, the quadrupole formalism is sufficient to derive 
the required information. It is known that two pointlike masses revolving around 
their common center of mass emit gravitational energy because of their time varying 
quadrupole moment. To describe their motion it is convenient to choose a coordinate 
system with the origin located at the center of mass of the system, the orbital plane 
coincident with the x-y plane and the x-axis oriented along the relative position 
vector, X = X* — Xp, when the planet is at the periastron. The vector X describes. 
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in general, an ellipse with semi-major axis a and eccentricity e, and its evolution is 
described by the parametric equations: 



a 



a{l — ecosu) 



2 arctan 



1-e 



1/2 



U 

tan — 
2 



(1) 
(2) 



Here p = |X|, a is the angle between X and the a;— axis, and u is the eccentric 
anomaly, related to time by the Kepler equation: 



Wkt = u — esinu, 
where iVk is the keplerian orbital frequency 



(3) 



_ 27r _ /GM\ 



1/2 



and M = + Mp is the total mass of the system. From the reduced quadrupole 
moment Qki, given by: 



(4) 



where fj, = M^,Mp/M is the reduced mass, it is straightforward to compute the 
amplitudes of the metric perturbation, projected on a sphere at distance r from the 
source: 



hij{t,r,l 



2G 



P,kPjl - \P^JPH ] Qklir) 



(5) 



r=t—r/c 



where Pij = (5* — riirij is the projector onto the plane orthogonal to the radial unit 
vector n. The explicit expressions of h0$ and he^ evaluated by this procedure 
are 



h00{u,r,0, 



c^ra(l — ecosu)^ 
X {(e^ - l)(l + cos2 6')(l -2cos2()!.) 

+ cos^ u(ecosu — 2) [(1 — e^)(cos^ </> + cos^ 0cos^ — cos^ 9) 
+ cos^ <f> + cos^ 6cos^ cf) — l\ + ecosu(cos^ (f> + cos^ 9cos^ — 1) 
+ 2\/l — sin^cos^(l + cos^ 9) sinu(ecos^ u — 2 cosu + e)}. 



(6) 



h0^{u,r,9,(j)) 



AG^Mn cos 9 
c^ra(l — ecosu)^ 

{sin (/)cos (p [(e^ — 2) cos^ u{e cosw — 2) — e cosu + 2(e^ — 1 
(2cos^ (j) — 1) vl — e^sinu(ecos^ u — 2cosw + e)}. 



(7) 
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If the orbit is circular the radiation is emitted at twice the orbital frequency 
(w*^^ = 2ujk). By Fourier-transforming the wave amplitudes, and by averaging 
over the solid angle, it is easy to show that if the orbit is eccentric, waves will be 
emitted at frequencies multiple of lu^ , and the number of equally spaced spectral 
lines will increase with the eccentricity. The average energy flux Fn relative to the 
n— th harmonic is given by: 



SttG 



(8) 



where (jig^g)'^ and (fi'J^)'^, are the square of the n— th Fourier component of the two 
independent polarizations, averaged over the solid angle 

= ^1'^^ \hee{nu;k,r,e,c^)\\ (h^^^'^f ^ h J '^^ \hg^{ncjk,r,e,cj,)\' . 

An estimate of the characteristic amplitude of the gravitational waves emitted 
by a planetary system can now be given by using the well known formulallj 

1/2 



K{nco,,r)^^±[{hl;;^f + {hi-Jr\ , (9) 

where the factor •y/2/3 takes into account the average over orientation. The afore- 
mentioned procedure, which is equivalent to that introduced by Peters and Mathewalj, 
has been applied to compute the characteristic wave amplitude he impinging on 
Earth, emitted by our set of EPS's and, for comparison, by the binary system PSR 
1913+16. Some results are shown in figure 1, where he is plotted as a function 
of the harmonic index n. It should be reminded that PSR 1913-1-16 is composed of 
two very compact stars with masses mi = 1.4411 M© and TO2 = 1-3874 Mq, 
revolving around their center of mass with an eccentric orbit (e 0.617139), semi- 
major axis a = 1.9490 • 10^^ cm, and keplerian frequency Vk = 3.583 • 10~^ Hz. 
The binary system is at a distance D = 5 kpc from Earth. In the upper panel of 
figure 1, we show two systems in which the planet moves in a nearly circular orbit 
around the central star, HD 283750 (e = 0.02) and r Boo (e — 0.018). As expected, 
the emission is concentrated at twice the keplerian frequency, and the amplitude is 
comparable to the maximum wave amplitude reached by PSR 1913-1-16, which is 
shown for comparison at the bottom of the figure. However, the frequency is about 
ten times lower than that of the binary pulsar. In the lower panel we show the 
characteristic emission of two EPS's with high eccentricity, 16 Cyg B (e = 0.634) 
and HD 89707 (e = 0.93). In this case the gravitational emission is spread over a 
larger set of frequencies, all multiple of i^fc, as for PSR 1913-1-16. 

In table 3 we tabulate the keplerian frequency, the maximum characteristic 
amplitude on Earth, and the corresponding frequency, for each EPS. In the last 
row the same data are listed for PSR 1913-1-16. For most systems the maximum 
emission frequency appears to be extremely low, in general smaller than 10^® Hz 
except the case of HD 283750, for which Vmax = 1-3 • 10~^ Hz. 
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HD 283750 
v-6 , 



10 12 14 16 18 20 

n 



hJO' 



23 




16 Cyg B 



10 12 14 16 



HD 89707 



Vi^=3.9 10'° Hz 



60 80 100 120 14 

n 



hJO'' 



0.8 - 

0.7 - 

0.6 - 

0.5 - 

0.4 - 



PSR 1913+16 



10 12 14 16 18 

n 



Figure 1: The characteristic amphtude computed from eqs. (||) is plotted versus the 
harmonic index n for four selected EPS's and for the binary system PSR 1913+16. 
For the two systems in the upper panel, HD 283750 and r Boo, the planet moves 
on a nearly circular orbit, with e = 0.02 and e = 0.018, respectively. In this case 
the emission is concentrated at twice the keplerian frequency, and the amplitude 
is comparable to the maximum wave amplitude emitted by the binary pulsar PSR 
1913+16. In the lower part of the figure he is plotted for two systems with high 
eccentricity: 16 CvK B (e = 0.634) and HD 89707 (e = 0.93), and for PSR 1913+16 
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Table 3: The maximum characteristic amphtude of the waves emitted by the se- 
lected set of EPS's and by the binary system PSR 1913+16, due to their orbital 
motion. 



Star 


[11 -) 










[Hz) 


HD 75289 


3.3 


10- 


-6 


2.1 


10- 


-25 


6.6 


10- 


-b 


51 Peg 


2.8 


10- 


-6 


4.0- 


10- 


-25 


5.5 


10- 


-6 


V And 


2.5 


10- 


-6 


7.9- 


10- 


-26 


5.0 


10- 


-6 


55 Cnc 


7.9 


10- 


-7 


3.8 


10- 


-25 


1.6 


10- 


-6 


p CrB 


2.9 


10- 


-7 


1.5- 


10- 


-25 


5.8 


10- 


-7 


HD 210277 


2.6 


10- 


-8 


1.9- 


10- 


-26 


7.9 


10- 


-8 


16 Cyg B 


1.4 


10- 


-8 


1.6 


10- 


-26 


5.8 


10- 


-8 


Gl 876 


1.9 


10- 


-7 


3.5- 


10- 


-25 


3.8 


10- 


-7 


47 Uma 


1.1 


10- 


-8 


5.3- 


10- 


-26 


2.1 


10- 


-8 


14 Her 


7.1 


10- 


-9 


2.7 


10- 


-26 


1.4 


10- 


-8 


Gl 86 


7.3 


10- 


-7 


1.6- 


10- 


-24 


1.5 


10- 


-6 


T Boo 


3.5 


10- 


-6 


6.6- 


10- 


-24 


7.0 


10- 


-6 


HD 168443 


2.0 


10- 


-7 


1.6 


10- 


-25 


6.0 


10- 


-7 


70 Vir 


9.9 


10- 


-8 


3.6- 


10- 


-25 


2.0 


10- 


-7 


HD 114762 


1.4 


10- 


-7 


2.0- 


10- 


-25 


2.7 


10- 


-7 


HD 110833 


4.3 


10- 


-8 


2.6 


10- 


-25 


2.1 


10- 


-7 


HD 112758 


1.1 


10- 


-7 


3.0- 


10- 


-24 


2.2 


10- 


-7 


HD 29587 


1.0 


10- 


-8 


2.3- 


10- 


-25 


2.0 


10- 


-8 


HD 283750 


6.5 


10- 


-6 


3.7 


10- 


-23 


1.3 


10- 


-5 


HD 89707 


3.9 


10- 


-8 


8.4- 


10- 


-25 


3.9 


10- 


-8 


HD 217580 


2.5 


10- 


-8 


8.5- 


10- 


-25 


7.6 


10- 


-8 


PSR 1913+16 


3.6 


10- 


-5 


1.0- 


10- 


-23 


1.4 


10- 


-4 
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4. Resonant excitation of the modes of a star by the orbiting planet 

We shall now consider another mechanism through which gravitational waves 
can be emitted by a system composed of a star and a planet. Since the masses 
of planets are much smaller than those of stars, in what follows we shall treat 
the planet as a pointlike mass, whereas the central star will be assumed to have a 
structure, and to possess a set of eigenmodes of oscillation, the quasi-normal modes, 
that are associated to the emission of gravitational radiation. We want to establish 
whether these modes can be excited by a planet. 

As mentioned in the introduction, perturbative calculations in the framework of 
general relativity made by KojimaEi have shown that, when a small mass moves in a 
circular orbit with orbital frequency Wfc around a compact star, the fundamental 
mode is excited if its frequency is — '2(-0k- In this case, a sharp resonance 

occurs, and the characteristic wave amplitude can be considerably larger than that 
evaluated by the quadrupole formula. It may be reminded that the quadrupole 
radiation by a planet in circular orbit around a star is also emitted at twice the 
orbital frequency, thus the condition of resonant excitation can also be formulated 
in terms of the quadrupole emission frequencies. 

Similar calculations have never been done for solar type stars; however, Ko- 
jima's results suggest that the resonant excitation of a mode may enhance the grav- 
itational emission also for non compact stars. For this reason, it is interesting to 
check whether the conditions of resonant excitation can be fulfilled in the planetary 
systems listed in table 1. We shall first verify whether the maximum quadrupole 
emission frequency of the planets of our set of EPS's, I'max (table 3, column 4) is 
close enough to any of the frequencies of the modes of the central star. 

The oscillation frequencies of a star can be computed if we make an assumption 
on its internal structure, i.e. on the equation of state prevailing in the interior. 
We shall consider, as an example, a very simple, polytropic model of star, with 
polytropic index n = 2. The oscillation frequencies ofjaewtonian polytropic stars 
are known to scale with the mean density of the starEII. In table 4 we tabulate 
the dimensionless eigenfrequencies of the modes, Vmodel {GM^, / R^)'^/'^ , for the 
chosen value of n. In order to explicitely compute the frequency of a given mode 
for a given star, the entries of table 4 have to be multiplied by the entries of table 
1, column 5. For instance, for the star HD 89707 the frequency of the /-mode is 
given by = 0.28 x 6.1 • 10^* = 1.7 ■ 10^* Hz. This number has to be compared 
whith Vmax = 2.2 • 10~^ Hz, given in table 3 for the same star, which is much 
smaller. This means that the planet cannot excite the /-mode of the central star. 
By repeating the same calculation for all systems and all modes, we find that it 
is unlikely that the stellar modes are excited in the EPS's we consider. This is 
because the angular velocities reached by the planets are too low, and consequently 
the frequencies of the radiation they emit are lower than those of the /-mode or of 
the lowest-order 5-modes of the star. Higher order g-modes could be excited, but 
the efficiency in producing gravitational radiation by this process would be too low. 

However, it should be reminded that according to recent theories on the evolution 
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Table 4: The dimensionless eigenfrequencies of the modes of a polytropic stars are 
tabulated for the polytropic index n — 2. 



mode i^mode/iGAh/Rjyr^ mode i^^ode/iGM^/Ri)^' 



PlO 


2.58 


/ 


0.28 


P9 


2.37 


91 


0.119 


PS 


2.15 


92 


0.087 


P7 


1.93 


93 


0.068 


P6 


1.70 


94 


0.056 


P5 


1.47 


95 


0.048 


Pi 


1.24 


96 


0.042 


P3 


1.01 


97 


0.037 


P2 


0.78 


98 


0.033 


Pi 


0.54 


99 


0.030 






910 


0.028 



of planetary systemsllil, there could exist planets moving on orbits even closer to 
the central star than those observed until now. In view of this possibility, it is 
interesting to investigate whether it is possible for a planet to approach a star at 
such a short distance that its angular velocity is high enough to excite the /-mode or 
the lowest-order g-modes, without being disrupted by tidal forces. In addition, we 
shall impose that the star does not accrete matter onto the planet. These conditions 
are equivalent to impose that neither the planet nor the star overflow their Roche 
lobe. We would like to stress that this limit takes into account only the gravitational 
interaction between the planet and the star. There may exist other processes that 
would prevent the planet to reach the innermost orbit allowed by the Roche-lobe 
analysis. However, as far as we know, the present knowledge on the formation of 
planets and on their possible migration toward the central star, still does not allow 
to firmly establish what is the minimum distance from a star at which a planet can 
safely sit. 

We shall assume, for simplicity, that the planet flies on a circular orbit with 
keplerian angular velocity ujk, and consequently emits radiation at the frequency 
ijjGW _ 2ujk. Let us indicate the dimensionless frequency tabulated in table 4 
multiplied by 27r as kmode ■ for instance, kg-^ = 0.750. The condition 
that a mode of the star is excited by the resonant interaction with the planet, 
i^mode = 2a;fe; therefore becomes 





1/2 

= 2 


'GiMp + A'hY 


1/2 


[ Rl \ 




[ a3 J 





(10) 



which can be written as 
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This equation gives the value the separation star-planet must have for a given mode 
to be excited. 

The further condition that the planet lies inside its Roche lobe, Rp < Rrl, is 
equivalent to the following constraint on its density: 



Pp > PRL, 

where orl = ^ — is the critical density. If we now introduce the dimensionless 

_ 3 ^^R.L 

quantity Rrl as given by 

Rrl = a Rrl, (12) 
i.e. we set to 1 the radius of the orbit, by the use of eq. ( pT| ) we find 



Prl 



which can be rewritten as 



.(1 



i-2 V ' AT 



Rrl 



(13) 



P* 4il + Mp/AQRj,L 

where is the mean density of the central star. Thus, a planet can excite a mode 
of the star corresponding to an assigned kmode, without overflowing its Roche 
lobe, if the ratio between its mean density and that of the central star exceeds the 
critical ratio (0). 

We have computed the dimensionless radius of the Roche lobe Rrl and the 
corresponding critical ratio (f4|), for assigned values of the ratio Mp/M^,, and 
the results are shown in table ||, which has to be read as follows. Suppose that the 
ratio between the mass of a planet and that of the central star is Mp/M^ = 10"'^. 
The frequency of the quadrupole radiation emitted by the planet will coincide with 
that of the first g-mode of the star, and the planet will not be disrupted, only if its 
density is higher that 1.49 p*. We have also checked whether the star overflows 
its Roche lobe and accretes matter onto the planet, and excluded from table ||the 
corresponding cases. It should be reminded that the ratio between the mean density 
of the planets of the solar system and that of the Sun is 3.9 for Mercury and the 
Earth, 3.7 for Venus, 2.8 for Mars, 0.9 for Jupiter, etc. A comparison of these values 
with the data of table || suggests that in principle, there can exist EPS's in which 
the first g-modes could be excited by a resonant process. For instance, a planet like 
the Earth could approach a polytropic star (n = 2) with the mass of the sun at a 
distance close enough to excite the mode gi without being disrupted by the tidal 
interaction, whereas a planet like Jupiter could only be at a distance good to excite 
the second 5-mode. 

By the Roche-lobe analysis we can also deduce another interesting information. 
Suppose that a planetary system made of a star with the mass of the Sun and a 
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Table 5: Ratios of the critical density of the planet pjiL to the central star density 
for different values of the mass ratio, Mp/Mi,, and different oscillation modes of 
a polytropic star with n — 2 (see text). 
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0.10 


0.10 


0.12 
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0.08 


0.08 


0.08 


0.08 


0.09 


0.10 



planet in circular orbit, is located at a fiducial distance D = 10 pc. We do not make 
any assumption on the internal structure of the star, but assign the values of the 
mass and of the mean density of the planet. In particular we consider four planets, 
with mass and density equal to that of Mercury, of the Earth, of Jupiter axid one 
with a mass equal to 13 times the mass of Jupiter and similar mean densityll3. We 
want to answer the following questions: 

• what is the minimum radius of the orbit? 

• what is the corresponding quadrupole emission frequency i/^^ = 2i/fc? 

• what is the corresponding characteristic amplitude on Earth? 

The answer is in table 6, where the required data are given for the four planets. 

From table 6 we see that planets like Jupiter or bigger could emit quadrupole 
radiation at a frequency in the bandwidth of space interferometers, and with an 
amplitude which could be even ten times bigger than that emitted by the binary 
pulsar PSR 1913+16. In addition, if the quadrupole radiation resonates with a 
mode of the star, the amount of emitted energy could be even larger. 

5. Concluding Remarks 

In this paper we have studied some of the EPS's discovered up to now, with the 
aim of characterizing their gravitational emission. 

As far as the quadrupole emission is concerned, we have shown that among those 
systems there is one, HD 283750, that emits a signal whith a maximum amplitude on 
Earth even higher than that of the binary pulsar PSR 1913+16, but at a frequency 
which is six times smaller. Moreover, since the orbit of the planet is nearly circular, 
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Table 6: We tabulate the minimum radius of the circular orbit, Qmin, the quadrupole 
emission frequency, z^*^^, and the characteristic amplitude of the corresponding 
wave emitted by four planets orbiting around a star with the mass of the Sun. We 
assume that the planets have mass and density equal to those of Mercury, Earth 
and Jupiter. The last planet has 13 times the mass of Jupiter and the same density. 
These data are obtained by imposing that the planets lie inside their Roche lobe 
(see text), and that the system is located at = 10 pc from Earth. 
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the radiation is almost entirely emitted at twice the keplerian frequency, whereas 
for systems with high eccentricity, as the binary pulsar, the radiation is emitted 
also at higher multiples of ujk- 

A further mechanism of gravitational emission is the resonant excitation of a 
mode of the star due to the tidal interaction with a planet. For compact stars, this 
mechanism has proved to be much more efficient than the quadrupole emission, thus 
it is interesting to check whether the conditions needed for the resonant excitation 
of the /- and g-modes can be fulfilled in any of the observed EPS's. Although there 
exist planets that have a very short orbital period (up to 1.79 days), and orbital 
frequencies close to 10~^ Hz, these values arc too low to induce a resonant 
excitation of the /-mode or the lowest ^r-modes of a solar type star. 

However, since many new EPS's will certainly be discovered in the near fu- 
ture, we wanted to understand whether more favourable conditions for the emission 
of gravitational waves may occur in these systems. In particular, since a higher 
emission frequency would be desirabile both for the excitation of the lowest order 
5-modes of the central star, and for a possible detection by space interferometers, 
we have investigated how close can a planet move on a circular orbit around a star 
without being tidally disrupted, and without accreting matter from the star. 

It should be stressed that the limits we establish by the Roche-lobe analysis 
do not take into account other processes which may destabilize the orbit of the 
planets, whose study is beyond the scopes of this paper. Having this caveat in 
mind, we have established that, in principle, there could exist systems in which 
the excitability conditions of the lowest order gr-modes could be fulfilled. In this 
case radiation would be emitted at frequencies of the order of ~ 10"^ Hz. The 
amplitude of the emitted gravitational signals depends on the ratio Mp/M^,, and 
for a planet like Jupiter or bigger, located at a distance of 10 pc, it would range 
between 10~^^ — 10~^^. The emitted radiation could be even larger, if the system is 
in a condition of resonant excitation of a mode of the star, and we plan to investigate 
this problem in detail in a subsequent paper. 



14 Gravitational waves emitted by extrasolar planetary systems 



Acknowledgements 

We would like to thank K. Kokkotas and L. Stella for their valuable suggestions. 
References 

1. A. Wolsczan and D. A. Frail, Nature, 355, 145 (1992). 

2. J. Schneider, Extrasolar Planets Catalog 
http : //www . obspm . f r/ encycl/catalog. 

3. R. A. Hulse and J. H. Taylor, Astrophys. J. 195, L51 (1975). 

4. J. H. Taylor and J. M. Weisberg, Astrophys. J. 253, 908 (1982). 

5. S. E. Thorsett and D. Chakrabarty, Astrophys. J. 512, 288 (1999). 

6. M. E. Alexander, MNRAS 227, 843 (1987). 

7. C. Terqucm, J. C. B. Papaloizou, R. P. Nelson and D. N. C. Lin, Astrophys. J. 502, 
802 (1998). 

8. K. D. Kokkotas and G. Schafer, MNRAS, 275, 301 (1995). 

9. Y. Kojima, Prog. Theor. Phys. 77, 297 (1987). 

10. http://obswww.unige.ch/ udry/planet/hd75289.htinl 

11. E. B. Ford, F. A. Rasio and A. Sills, Astrophys. J. 5 14, 411 (1999). 

12. G. M. Marcy, P. Butler, S. Vogt, D. Fischer and J. Lissauer, Astrophys. J. 505, L147 
(1998). 

13. M. Mayor, D. Quoloz, S. Udry and J. L. Halbwachs, in Bioastronomy 96, Astronom- 
ical and Biochemical Origins and the Search for Life in the Universe, ed. C. B. 
Cosmovici, S. Bowyer and D. Werthimer (1997). 

14. A. Duquennoy and M. Mayor, Astron. Astrophys. 248, 485 (1991). 

15. K. S. Thorne, in 300 Years of Gravitation, ed. S. Hawking and W. Israel (Cambridge 
University Press, Cambridge, 1987). 

16. P. C. Peters and J. Mathews, Phys. Rev. 131, 435 (1963). 

17. K. P. Cox, Theory of stellar pulsations (Princeton: Princeton University Press, 1980). 

18. D. Saumon, W.B. Hubbard, A. Burrows, T. Guillot, J.I. Lunine and G. Chabrier, 
Astrophys. J. 460, 993 (1996). 



